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STOCHASTIC  ESTIMATION  IN  COMBINED  ARMS 
LANCHESTER  MODELING  OF  WARFARE 
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^  Deterministic  Lanchester  models  of  wajcfare  with  time  dependent  attrition 
coefficients  has  received  extensive  study'll] V  There  does  not  appear  to  be 
any  analytical  effort  to  incorporate  observations  of  the  battlefield  into 
Lanchester  models.  This  is  very  surprising  when  one  considers  the  vast  amount 
of  effort  going  on  in  developing  techniques  for  observing  events  on  the 
battlefield  and  integrating  intelligence  in  a ^fusion'^enter.  In  addition, 
no  analytical  work  appears  to  have  incorporated  uncertainties  in  the  model 
itself.  The  objective  of  this  paper  is  to  develop  a  general  methodology, 
using  stochastic  estimation  theory  as  practiced  in  modern  control  theory,  for 
predicting  conditions  on  the  battlefield  from  measurements  or  observations  of 
the  battle  area  that  are  themselves  uncertain  and  where  the  parameters  of  the 
mathematical  model  are  not  known  precisely.  Examples  are  also  provided  to 
show  how  the  general  theory  can  be  applied  to  two  specific  problems. 

/The  fundamental  Lanchest^equ^ions  for  homogeneous  forces  is  given  by 
dx 

—  =  -  ay 

dt 

(1) 

dy 

—  =  -  gx 
dt 

where  x  and  y  are  the  force  strength  of  each  side  and  a  is  the  attrition  of  x 
force  due  to  the  weapon  effectiveness  of  the  y  force  and  similarly  p  is  the 
attrition  coefficient  of  the  y  force  due  to  the  weapon  effectiveness  of  the 
X  force.  Taking  into  consideration  that  the  values  a  and  p  can  change  with 
time  and  also  are  not  known  exactly,  we  can  incorporate  the  uncertainties  in 
these  parameters  by  adding  to  our  model  a  gaussian  white  noise  term  in  each 
equation  so  that  equation  (1)  is  modified  as 

dx 

—  =  -  ay  +  wi  (t) 
dt 

(2) 

dy 

—  =  -  px  +  W2  (t) 
dt 

where  wj  and  W2  are  uncorrelated  gaussian  white  noise. 

Equation  (2)  can  be  written  in  matrix  form 

jC)  ■("  •)(:) 

or  with  =  (x,  y)T  vector 


0 


matrix 


A  = 


-a 


§ 


-e 


w  =  (w^,  W2)'*’ 


vector 


we  have 

dx 

dt 


A  X  +  w 


(3) 


where  w(t)  is  gaussian  white  noise  such  that 
E[w(t)]  =  0  for  all  t 

cov  [w(t) ,  w(t)]  =  ^5  (t-x)  (3a) 

i  .e. ,  w  ~  N  (0, 
with  ^  sO 

Equation  (3)  and  (3a)  represents  the  vector  form  of  the  homogeneous  Lanchester 
model  with  "plant"  noise  w 

Observation  over  the  battlefield  can  be  similarly  characterized  as 

zi  =  Cl  X  +  VI 

Z2  =  C2  y  +  V2 

where  zi  is  the  measurement  of  x  corrupted  by  measurement  noise  vj,  and  Z2  is 
the  measur^ent  of  v  corrupted  by  measurement  noise  V2.  If  we  let  z-  (z^,  Z2)', 
Ji  ~  (Xf  y)  and  *  (vj,  V2)^  we  have 


(4) 


Cl  0 
0  C2 

and  ^(t)  ~  N  (0,  ^  (t))  gaussion  white  noise 
such  that 

E  I  ^  (t)  I  =  0  for  all  t 
cov  I  ^  (t),  1  (t)  I  =  R  5  (t-x) 

R(t)  »  RT(t)  >0 


(4a) 


2 


We  also  require  R  to  be  non-singular  i.e.,  R"1  exists. 

Thus  the  fundamental  set  of  equations  which  describe  the  observations 
over  the  battlefield  for  “stochastic"  homogeneous  Lanchester  models  of  warfare 
is  given  by 

dx(t) 

(System  model): - *Aii(t)  +  w(t)  ) 

dt  I  (5) 


(Measurement  model):  ^(t)  =  x.  (t)  +  (t)  | 

where  the  initial  state  2c(to)  is  modeled  as  a  random  vector  since  the  initial 
conditions  are  not  precisely  known.  To  characterize  the  initial  state  x  (tg) 
we  need 

E  {x  (to)  I  =  io 

cov  [  X  (to),  X.  (to)]  =  £0  initial  state  pvariance  (6) 

P  a  pT 

Io  i-o  ^ 

The  "plant"  noise  w(t)  is  used  to  model  external  disturbances  and  modeling 
errors  in  A.  Similarly,  the  measurement  noise  y  (t)  is  used  for  modeling 
error  in  C. 

It  is  also  assumed  that  ^((to),  w(t),  and  y(t)  are  independent  for  all  to,  t, 
T,  i.e., 

cov  lx  (to),  w  (t)]  =  0  V  to,  t 

cov  [x  (to),  y  (t)]  =  0  V  to,  t  (7) 

cov  [w  (t),  y  (t)]  =  0  V  to,  t 

; 

The  flow  diagram  for  equation  (5)  is  shown  in  Figure  1. 

Thus,  the  stochastic  estimation  problem  is  to  find  an  estimate  £  for  the 
stochastic  homooeneous  Lanchester  equation  (5)  subject  to  equations  pa), 
(4a),  (6),  and  (7)  such  that,  for  a  given  set  of  past  measurements  i(l), 
z^(2),  ...  ,  ^(t),  the  estimate  £(t)  minimizes 

J  =  E  I  ||x(t)  -  x(t)l|2  I  z(l),  z(2),  ...  ,  z(t)|  (7a) 

in  the  least  squares  sense. 

It  is  now  well  known  that  the  Kalman-Bucy  filter  [2,  3,  4]  provides  such  an 
estimate  and  that  x  satisfies 


(9) 


where  the  filter  gain  matrix  H  is  given  by 
H  =  P(t)  CT  R-1  (t) 

with  the  error  covariance  matrix  £(t)  satisfies  the  Riccati  equation 


dP(t) 

dt 


A  P(t)  +  P(t)AT  +  -  P(t)  CT  R-1  C  P(t) 


(10) 


with  P(to)  =  £0 

The  Kalman-Bucy  filter  can  be  implemented  as  shown  in  Figure  2. 

The  state  estimate  ^  provides  a  battlefield  commander  with  a  predicted  estimate 
of  the  future  conditions  of  his  force  strength  and  that  of  his  opponents  based 
on  available  information  up  to  the  present  time  for  a  homogeneous  force. 


The  heterogeneous  or  combined  arms  Lanchester  model  of  warfare  can  be 
simply  achieved  by  expanding  the  definition  of  x  and  y  to  a  force  structure 
made  up  of  several  distinct  weapon  systems  x-j,  yj  for  i  =  l,2,  ...,n,j  = 
1,  2,  n  such  that 


X 


and 


w 

X 

z 


X2,  Xp, 

( Wj^  ,  W2  t  •  •  •  *  Wp  , 

(vj,  V2,  ....  Vp, 
(^1 ,  Z2 ,  . • • ,  Zp  , 


yi,  yz . ym)''^ 

'^n+1  ^n+m^ 

''n+1  *  ••*»  ''n+m^^ 
^n+1 »  •••*  ^n+m^ 


n  +  m  vector 

n  +  m  vector 
n  +  m  vector 
n  +  m  vector 


and  where  the  plant  matrix  F  takes  on  the  form 


F  = 


(n  +  m)  X  (n  +  m)  matrix 


where  A,  £  are  now  matrices  which  represent  the  attrition  coefficients  asso¬ 
ciated  with  weapon  yj  and  x-j  for  j  =  1,  2,  m  and  i  =  1,  2,  ...  n. 


KALMAN -BUCY  Filter 


Thus  for  the  heterogenous  case  equation  (5)  is  transformed  into 


dx 

—  =  F_  X 

dt  (11 

z{t)  =  C  X  +  V 

Naturally,  C  is  similarly  defined.  The  implementation  is  similar  to 
that  shown  in  Figure  1.  All  the  constraints  imposed  on  w  ^nd  ^  are  the  same 
as  in  the  homogeneous  case.  The  optimal  estimate,  generated  by  the  Kalman- 
Bucy  filter  for  the  heterogeneous  systems,  is  also  satisfied  for  equation 
(8)  through  (10).  The  Kalman-Bucy  filter  approach  can  also  be  applied  to 
the  situation  in  which  F(t)  and  C(t)  are  time  dependent. 

Discrete  System 

We  can  transform  the  homogeneous  system  (equation  (5))  to  discrete  true 
model  if  we  use 


d2<(t) 

dt 


x(tk)  -  x(tK.i) 

1  in  —  .1.1 

At->0  At 


so  that 

x(t^/)  -  x(tK.i)  =  X  (t^.i)  +  At  w  (t^.i) 
z(tk)  =  C  X  (tk.i)  +  X  (tk-i) 


which  reduces  to 

ii(tk)  =  *K-1  2i(tK-l)  +  At  w  (tK-i) 
l(tk)  =  i  ii  (k)  +  1  (tk) 
where  -\ 


(12 


♦k-1  ®  1 


WK-1  =  Atw(tK.i) 
w(tk)  ~  N  (0,  3k) 

5k  =  S(t)/At 
v(tk)  ~  N  (0, 


(13 


^  =  R(t)/At  ^ 

with  initial  state  2i(o)  gaussian  and  2i(o)»  w(tk-i),  1  (^k-l)  ^re  independent 
for  all  t ,  T  . 


The  opti3'il  filtering  problem  is  to  determine  x(tK_i)  given  a  sequence 

of  past  measurements  ^(1),  ^(2) .  z{tK_i)  where  the  optimal  estimate  is 

defined  as  the  conditional  mean 

(known)  x  (t^.i  1  tK.i)  =  E|  x  (tK-i)  1  z(l),  z{2) . z(tK-i) |  (14) 

where  p(j(  (t^^.i)  ]  ^(1),  Z2(2),  2^(tK_i)),  the  conditional  density  function, 

is  gaussian  and  known.  The  conditional  covariance 

(known)  P(tK.i  1  tK-i)  =  cov  Cx(tK-l) ,  x(tK-i)  1  2(1),  z(2) . z(tK.i)]  (15) 

We  are  interested  in  predicting  x^(tK)  given  observations  ^(1),  z[Z) , 

...,  l(tK>i)  before  measurement  z^(tK)  is  made  and  then  to  improve  that  estimate 
x(tK)  based  on  expanded  information  z(l),  z(2),  ...,  z(tK_i),  z(tK)  after  we 
have  made  a  measurement  or  observation  at  tK  i.e.,  i(t|(). 

The  discrete  Kalman  filter  can  be  summarized  as  follows: 

State  Dynamics  Model: 

x(tK)  =  1  K-1  X  (tK-l)  +  w  (tK-l) 

w(tK)  ~  N(0,  3k) 

Measurement  Model 

l(tK)  =  C  X  (tK)  +  V  (tK) 

v(tK)  ~  N  (0,  Rk) 

Off-Line  Calculations 
Initialization  (t=0): 

fo  =  E  [(x(0)  -  x(0))  (x(0)  -  x(0))T] 

Predict  Cycle: 

P.  (Tk  1  ^K.j)  =  1  (^K-l)  -(^K-l)  l^(^K-l)  S.(^k-l) 

Update  Cycle: 

-^^kI^k)  ”  -  ^^kI^k-i)  ”  -  (^k!^k-i)  (^k)  • 

[£(tK)  P.  (TkITk-i)  £^(tK)  ^  £(^k)  — ^^kI^k-i) 

Filter  Gain  Matrix: 

Hdn)  =  P(tK|t^)  cT(tn) 


The  structure  of  discrete  time  Kalman  filter  is  shown  in  Figure  4 


KALMAN -BUCY  Filtef 


•  X  •>»  M  •X  •>> 


Several  examples  will  now  be  considered  to  demonstrate  the  method  of  computa¬ 
tion,  details  of  which  are  presented  in  the  Appendix. 

Example  1:  Square  law  with  noise  in  only  one  of  the  state  variables. 

=  Ox  -  ay 

=  Bx  +  Oy  +  1  w  (noise  in  y  only) 

=  y  +  V  (observation  of  y  only) 


H  =  (0  1  ) 

E  (w(t'))  =  0  E  (v(t))  =  0 


E  (w(t)  w(t'))  =  q  6  (t-t')  E  (v(t)v(t'))  =  r  6(t-t') 


5  =  q  R  =  r 


E  (x  (0))  *  0;  E  (y  (0))  =  0 


Pl2  (0) 
P22  (0) 


P(0) 


P  =  A  P  +  P  aT  +  G  S  G T  -  P  hTr-1  h  P 


P21  *  P12 


(Riccati  Equation) 


12 


•<x|  o  •/x|  <x 


Pll  =  -2aP21  -  (  - 
\r/ 

)P21^ 

Pn  (0) 

given 

P12  =  -aP22  -BPii 

-^-^21  P22 

P12 

given 

P22  ■  -26Pi2  ■*■  P  ■ 

0-^ 

P22 

given 

Steady  State: 

^11  =  ^12  “  ^22  =  ° 
Pll  =  0  =>  P21  =  0 


.  1 

^12  (<*■•■“  ^21^  ^22  “  ^^11 

r 

•  ^  ? 

P22  *  0  =>  0  =  -26  P12  ^  -  ^22^ 

r 

with  Pj2  =  P21  =  0  =>  P22  “ 

P„  = 

=  A  X  +  G  w  +  P  hTr-I  (z  -  H  x)  X  (0)  =  0 

=  (A-K  H)x+Kz  K=P  hTr-1 


H  =  (0  1) 


A 

X 


A 

=  -ay 


P21  (t)  (Z  - 


y) 


where  x  (0)  =  0 


X 


at  steady  state 


5  J-  \  P22(t)  (z  -  5)  y  (0)  =  0 

Yq\l/2 

=  -b5  (z  -  P)  at  steady  state 

Summarizing  the  results  for  the  case  where  only  one  state  variable  (y) 
is  corrupted  by  noise  in  the  systems  model  and  measurement,  the  general 
solution  satisfies  the  equations 

"x  =  -oy  +  (P2i(t)/r)  (z  -  y)  x(0)  =  0 

y  =  -  Bx  +  (P22(t)/'")  (z  -  y)  y(0)  =  0 

while  the  steady  state  equations  take  the  form 

X  =  -ay  x(0)  =  0 


•  q  1^2 

y  *  -Bx  +  -  (z  -  y)  y(0)  =  0 

r 


In  the  next  example  (example  2)  we  consider  the  case  where  the  "plant"  ntise 
and  measurement  noise  corrupts  both  state  variables  as  follows. 

Example  2: 


X  =  Ox  -  ay  +  1  wj. 
y  =  -Bx  +  Oy  +  1  W2. 
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Pll  = 


/  ^ 
/  R?? 


Steady  state 


Steady  state 


Pl2  =  P21  0 

The  predicted  state  updated  by  measurements  at  an  earlier  time  now  is  given  by 


X  =  A  5  +  6  w  +  P  hT  R-1  (z  -  H  x) 


x(0)  =  0 


x=  (A-K)  X  +  K  z 


1  h\\  ^12  \ ^ 


R|  \  C21  522  /\  Z2  -  y 


X  (t)  =  -a  y  (t)  +  -  I  Cii  (t)  (zi  -  x)  +  Cl2(t)  (Z2  -  y) | 

1R| 


y  (t)  =  -B  X  (t)  +  -  I  ?2l(t)  (^1  -  x)  +  C22(t)  (^2  "  y)  [ 

|R| 

where 

5ii(t)  =  Pii(t)  R22  -  Pi2(t)  R21 

Cl2(t)  =  Pi2(t)  Ril  -  Pii(t)  Ri2 

52l{t)  =  P2l(t)  R22  -  P22(t)  R2I 

C22(t)  =  P22(t)  Pll  -  P2l(t)  Ri2 


Accordingly  we  have  obtained  either  the  differential  equation  for  which  the 
predicted  estimate  x  needs  to  satisfy  or  the  integral  solution  for  this 
estimate  in  terms  of  its  transition  matrix  ♦  (t,  to). 

In  summary,  there  currently  exists  engineering  techniques  from  the 
field  of  modern  control  theory  that  offers  the  commander  in  the  field  the 
opportunity  of  predicting  or  forecasting  conditions  of  the  force  strength  on 
the  battlefield  given  prior  information  or  observation  on  status  of  battle. 

The  approach  can  handle  any  number  of  different  weapon  systems.  The  methodology 
described  in  this  report  will  be  applied  to  recently  acquired  data  on  tank 
duels  and  will  be  reported  in  a  separate  report. 
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APPENDIX 


Detailed  analysis  and  derivations  will  now  be  presented  to  demonstrate  the 
method  of  computation. 

Example  1:  Square  law  with  noise  in  only  one  of  the  state  variables. 

=  Ox  -  ay 

=  3x  +  Oy  +  1  w  (noise  in  y  only) 

=  y  +  V  (observation  of  y  only) 


H=  (0  1  ) 

E  (w(tO)  =  0  E  (v(t))  =  0 

E  (w(t)  w(t'))  =  q  6  (t-t')  E  (v(t)v(t'))  =  r  6(t-t') 


23 


p  = 


(0)) 

=  0; 

z 

'  X 

L  P  + 

p  aT 

Pll 

P12 

P2I 

P22 

E  (y  (0))  =  0 


Pll  (0)  Pi2  (0) 


21 


(0)  P22  (0) 


P21  =  P12 


0  -c 


A  P  = 


Pll  P12 


-6  0  /  \  P21  P22 

-aP2i 


-aP21  -aP22 

-6P1I  -PP12 


(A  P)T  =  P  aT 


-“P22 


-CPll 

-BP12 


G  0  gT  =  q  /  °  Vo  1  j 


0 

0 


0 

1 


H  P  =  (0  1)  /Pll  P12 

P2I  P22 


(P2I.  P22) 


1  /  P 


(P21,  P22)  =  - 


21 


r  I  P21F 


’12  '  ‘“^22  '^^11  •(  ■  r21  ^22 


P22  =  -2aPi2  +  Q  -  -  P22 


=  P(0) 

(Riccati  Equation) 


Pll  (0)  given 
P12  (0)  given 


P22  (0)  given 


steady  State: 


^11  "  ^12  '  ^22  =  ° 


Pll  =  0  =>  P21  =  n 


Pj2  =  0=>0  =  -  (a+  -  P21)  P22  “  fiPii 

r 


^22  ~  0  0  =  “2B  P22  +  Q  ~  ~  ^22^ 

r 


with  Pj2  "  ^21  “  ^  ^22  " 


^11  "  ■  *  ^22 


-hrq)l/2 


x  =  Ax  +  Gw+P  hVI  (2  - 


X  (C)  «  0 


X  -  ( A-K  H)  £  +  K  2^ 


K  =  P  hTr-I 


'■(;—(.:  ■; 


H  =  (0  1) 


•  1  ^21  -  y) 


where  x  (0)  =  0 


at  steady  state 


..,X  *  - 


P22(t)  (2  -  y) 


(z  -  y) 


y  (0)  =  0 


at  steady  state 


Sunmarizing  the  results  for  the  case  where  only  one  state  variable  (y) 
is  corrupted  by  noise  in  the  systems  model  and  measurement,  the  general 
solution  satisfies  the  equations 


X  =  -ay  +  (P2i(t)/r)  (z  -  y) 


x(0)  =  0 


y  =  -  ex  +  (P22(t)/r)  (z  -  y) 
while  the  steady  state  equations  take  the  form 


y(0)  =  0 


A  A 

X  =  -ay 


^  =  -3x  +  (  —  )(z  -  y) 


x(0)  =  0 


y(o)  =  0 


In  the  next  example  (example  2)  we  consider  the  case  where  the  "plant"  noise 
and  measurement  noise  corrupts  both  state  variables  as  follows. 

Example  2: 


X  =  Ox  -  ay  +  1  w^, 


y  =  -6x  +  Oy  +  1  w;>, 


zi  =  X  +  VI 

Z2  =  y  +  V2 


0  V  >  0 

-6  0/  Vo  1 


=  A  X  +  I  w 


=  I  X  +  V 


E  (x(0))  =  0 


H  =  X;  V  = 


i' 

J 


E(w)  =  0;  E  (w(t)  w{t')T)  =  a(t)  6(t-t') 


,  ,  /  Qll  Q12  \ 

Q(t)  =  I  I  assumed  constant 

V  Q21  Q21  / 


E(v)  =  0;  E  (v(t)  v(t')''‘)  =  i(t)  6  (t-r);  R(t) 


R-1  -  i.  ( 


1R1\-R21  Rii 

R  1s  assumed  constant 

IRl  ♦  0 


’’‘2 

\P21  P22 

G  ^  gT  =  ^ 

H  P  =  P 


IRI  =  Rii  R22  -  R21  Ri2 


R-1  P  = 


R22  -R12 

"R21  Rll 


1  /  Pii  R22  -  Ri2  P2I 

1R|  \  P21  Rll  -  R21  Pii 


R22  P12  -  R12  P22 
Rll  P22  -  R21  P12 


■  •  -  f.  W  r  W  7*  V  •.  V 


p  R-lp  =  - 

IRI 


1  /  Pll  Pl2  W  ^11  ^^2  -  Pl2  P21 
PI  I  P21  P22  /  ^21  Pll  -  P21  Pll 


P22  P12 

Pll  P22 


P12  P22 
P2I  P12 


1  /an  ai2 

PIU2I  ^22 


where 


an  =  Pll  (Pll  R22  -  P12  P21)  +  P12  (P2I  Pll 

ai2  =  Pll  (R22  P12  -  P12  P22)  +  P12  (Pll  P22 

321  =  P21  (Pll  P22  -  P12  P21)  +  P22  (P2I  Pll 

322  =  P12  (P22  P12  -  P12  P22)  +  P22  (Pll  P22 


R2I  Pll) 
P2I  P12) 
P2I  Pll) 
P2I  P12) 


for  Ri2  =  R21  and  P21  =  P12  it  follows  that  ai2  =  a2i. 
K  =  P  hT  R-1 


|P|  \P21  P22/ 

/Pll  P22  -  P12  P2I 


1R| 

\P21  R22 

CM 

CM 

0. 

1 

1 

' 

512  \ 

\521 

522/ 

where 

^11  =  Pll 

P22  -  P12 

P2I 

^12  =  P12 

Pll  -  Pll 

P12 

^21  =  P2I 

P22  -  P22 

P2I 

^22  =  P22 

Pll  -  P2I 

P12 

Pll/ 

Pll  P12  -  P12  Pll 

Pll  P22  -  P21  P12 


••-  W* 


-  :  v  v  .  'V 


kT  =  R-1  P 


P  =  AP  +  PAT'  +  fifiG'>’-KHP 
P(0)  =  Po 

KHP=P1R-1  IP 
=  P  R-1  P 

The  Riccati  equation  is  given  by 


1 


P  =  /  -2a  P21  Qll  "  ~  3ll 

I''  I 


-gPjj  -aP22  +Q12 


1 


-CPll  -aP22  +Q21  -  321  -5^12  -“^22  +Q2I 

1R| 


which  for  steady  state  £  =  0  and  P12 
Evaluating  £  at  steady  state 


1  /dll  «>12 

R-1  £  =  —  ( 

|R|  Vd2i  d22 


where 


^11  =  R22  Pll  -  R12  P2I 
di2  =  R22  P12  -  P12  P22 
<^21  =  Rii  P21  -  R21  Pll 
d22  =  P22  Pll  -  R2I  P12 


p  R-1  p  =  — 


^/Pll  Pl2\/<lll  ^512 

|P|\P21  P22/\^21 


1  /  Pll  dll  <i\Z  P12 

|R|  \P2i  dll  +  d2i  P22 


Pll  di2  +  d22  P12 


P21  di2  +  d22  P22  f 


29 


with  ?i2  — 0  dll  “  ^22  ^11 


di2  =  -Ri2  P22 


at  Steady  state: 


1  /  R22  ^11 

PR-lp=  22  11 

|Rl  I  -R21  Pii  P22 

So  that  at  steady  state 


(1)  0  =  -2  a  P21  +  Qii  -  ail 

IRI 


(2)  0  =  -e  Pii  -a  P22  +Q12  -  7“-  ai2 

|Rl 


(3)  0  =  -6  Pii  -a  P22  +  Q2I  -  7—  321 

IR! 


(4)  0  =  -2$  P12  +  Q22  ■  “  ^22 


Steady  state:  P12  >  0 

hi  — ^  ’’11^  ^22 

312  — >  -R12  P22  Pll 

321  - ■>  -R2I  Pll  P22 

322  — ^  '^ll  ^22^ 

1  o 

(la)  0  =  Q11  -  —  R22  Pn 
11 


(2a)  0  *  -2Pii  -aP22  +  Ql2  *  -  Pll  P22 

IRI 


d2i  =  -R21  Pll 
d22  =  Rll  P22 


■Ri2  P22  Pll 
Rll  *^22^ 


.  •-  v;.  -v . 


i 


i 

i 


R 


(3a) 

0  =  -pPii  -aP22 

12 

+  Q21  ~  Pii  P22 

|R! 

(4a) 

1 

0  =  Q22  ■ 

|RI 

Rll  P22^ 

Thus 

p  2_  “ 

Rll  “ 

R22 

IRI 

(la')  Pii 

=  /  On 

/  R22 

IRI 

steady  state 

(4a')  P22 

=  /  Q22 

/  Rll 

|R| 

Steady  state 

P12  =  P21  — > 

0 

The  predicted  state  updated  by  measurements  at  an  earlier  time  now  is  given  by 

• 

X  =  A  2  +  w  +  P  hT  R-1  (z  -  H  £)  x(0)  =  n 

or 


X  (t) 


•ay  (t)  + 


1R| 


]4ll  (t)  (zi  -  x)  +  4i2(t)  (Z2  - 


31 


-p  X  (t)  +  —  |C2l(t)  (zi  -  x)  +  C22(t)  (Z2  -  y)| 

IRI 


5ll(t)  =  °ll(t)  R22  -  Pl2(t)  R2I 
^12(t)  =  Pl2(^)  ^11  ■  ^12 
C2l(t)  =  P2l(t)  R22  -  P22(t)  P2I 
C22(t)  =  P22(t)  Rll  -  P2l(t)  Ri2 


■5ll(t) 


5l2(t) 


52l(t) 


2(t)  * 


-6  + 


^2l(t:) 


=  f2(t)  Ji  +  1  z_ 


cii(t) 


52l(t) 


7  - 


^12(t) 


522(t) 


y  + 


5l2(t) 


52l(t) 


C22(t) 


^zzit) 


Cl2(t) 


C22(<^) 


=  A(t)  -  K(t) 


=  K(t) 


^(t)  =  ♦  (t,  to)  x(to)  +  J  i  (t,  t)  K(t)  z(t)  dT 


where  *  (t,  t©)  satisfies  the  homogeneous  equation 


i  (t.to)  =  £2  (t)  *  (t,to) 

Accordingly  we  have  obtained  either  the  differential  equation  for  which  the 
predicted  estimate  ^  needs  to  satisfy  or  the  integral  solution  for  this 
estimate  in  terms  of  its  transition  matrix  «  (t,  to). 
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